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F^ i In this paper we present a self-contained macroscopic description of 

Ch ' diffusive systems interacting with boundary reservoirs and under the 

O ■ action of external fields. The approach is based on simple postu- 

lates which are suggested by a wide class of microscopic stochastic 
models where they are satisfied. The description however does not 
^ ' refer in any way to an underlying microscopic dynamics: the only in- 

C"^^ , put required are transport coefficients as functions of thermodynamic 

variables, which are experimentally accessible. The basic postulates 
■t:^ ■ are local equilibrium which allows a hydrodynamic description of the 

f-~^ . evolution, the Einstein relation among the transport coefficients, and 

C3 I a variational principle defining the out of equilibrium free energy. As- 

sociated to the variational principle there is a Hamilton-Jacobi equa- 
tion satisfied by the free energy, very useful for concrete calculations. 
Correlations over a macroscopic scale are, in our scheme, a generic 
^ ' property of nonequilibrium states. Correlation functions of any order 

j^ ■ can be calculated from the free energy functional which is generically 

a non local functional of thermodynamic variables. Special attention 
is given to the notion of equilibrium state from the standpoint of 
nonequilibrium. 
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1 Introduction 

The main purpose of the present paper is to describe a "physical theory" 
for a certain class of thermodynamic systems out of equilibrium, which is 
founded on and supported by the analysis of a large family of stochastic 
microscopic models. Out of equilibrium the variety of phenomena one can 
conceive makes it difficult to define general classes of phenomena for which 
a unified study is possible. Furthermore the details of the microscopic dy- 
namics play a far greater role than in equilibrium. Since the first attempts 
to construct a non equilibrium thermodynamics, a guiding idea has been 
that of local equilibrium, which means that locally on the macroscopic scale 
it is possible to define thermodynamic variables like density, temperature, 
chemical potentials... which vary smoothly on the same scale. Microscopi- 
cally this implies that the system reaches local equilibrium in a time which is 
short compared to the times typical of macroscopic evolutions, as described 
for example by hydrodynamic equations. There are important cases however 
where local equilibrium apparently fails like aging phenomena in disordered 
systems due to insufficient ergodicity. These will not be considered in this 
paper. Also the case in which magnetic fields play a role is not covered by 
our analysis. 

The simplest nonequilibrium states one can imagine are stationary states 
of systems in contact with different reservoirs and/or under the action of ex- 
ternal (electric) fields. In such cases, contrary to equilibrium, there are cur- 
rents (electrical, heat, matter of various chemical constitutions ...) through 
the system whose macroscopic behavior is encoded in transport coefficients 
like the diffusion coefficient, the conductivity or the mobility. 

The ideal would be to approach the study of these states starting from 
a microscopic dynamics of molecules interacting with realistic forces and 
evolving with Newtonian dynamics. This is beyond the reach of present 
day mathematical tools and much simpler models have to be adopted in the 
reasonable hope that some essential features are adequately captured. In the 
last decades stochastic models of interacting particle systems have provided 
a very useful laboratory for studying properties of stationary nonequilibrium 
states. From the study of these models has emerged a macroscopic theory for 
nonequilibrium diffusive systems which can be used as a phenomenological 
theory. 

A basic issue is the definition of nonequilibrium thermodynamic functions. 
For stochastic lattice gases a natural solution to this problem has been given 
via a theory of dynamic large deviations (deviations from hydrodynamic 
trajectories) and an associated variational principle leading to a definition of 
the free energy in terms of transport coefficients. 



One of the main differences between equilibrium and nonequilibrium sys- 
tems, is tliat out of equilibrium the free energy is, in general, a non local 
functional thus implying the existence of correlations at the macroscopic 
scale. These correlations have been observed experimentally [11] and appear 
to be a generic consequence of our variational principle which can be refor- 
mulated as a time independent Hamilton- Jacobi equation for the free energy. 
This is a functional derivative equation whose independent arguments are 
the local thermodynamic variables and which requires as input the transport 
coefficients. The presence of long range correlations for out of equilibrium 
states had been derived within kinetic theory, see [11] and reference therein, 
as well as within fluctuating hydrodynamics, see e.g. [25]. 

We believe that the theory we propose is a substantial improvement 
with respect to the theory developed long ago by Onsager [22] and then by 
Onsager-Machlup [23] which applies to states close to equilibrium, namely for 
linear evolution equations, and does not really include the effect of nontrivial 
boundary reservoirs. In principle, the theory we suggest should be applicable 
to real systems, i.e. with nonlinear evolution equations and arbitrary bound- 
ary conditions, where the diffusion is the dominant dynamical mechanism. 
We emphasize however that we assume a linear response with respect to the 
external applied field. Of course, the Onsager theory is recovered as ffist 
order approximation. 

The basic principle of the theory here proposed is a variational principle 
for the nonequilibrium thermodynamic functionals. This principle has the 
following content. Take the stationary state as the reference state and con- 
sider a trajectory leading the system to a new state. This trajectory can be 
realized by imposing a suitable additional external field (in addition to the 
one already acting on the system). Then compute the work done by this 
extra field and minimize it over all possible trajectories leading to the new 
state. This minimal work is identified with the variation of the free energy 
between the reference and the final state. As well known in thermodynamics, 
in the case of equilibrium states this definition agrees with the standard one. 

Our treatment is based on an approach developed by the authors in the 
analysis of fluctuations in stochastic lattice gases [1-7] . For recent overviews 
on nonequilibrium phenomena see also [13, 17, 21] and [24] for a different 
approach which includes a discussion on long range correlations. 

The outline of the paper is the following. 

In Section [2] we introduce the basic principles for a self-contained macro- 
scopic description of driven diffusive systems out of equilibrium. These sys- 
tems are described by local thermodynamic variables which evolve according 
to (nonhnear) diffusion equations with boundary conditions and/or external 



fields which force the system out of equilibrium. We emphasize that such 
evolution equations have the structure of conservation laws. The possibility 
of such description is of course based on the assumption of local equihbrium. 
We then formulate mathematically the above mentioned variational principle 
for the nonequilibrium free energy and characterize the optimal trajectory. 
As we already remarked, our approach gives the usual equilibrium free energy 
when we consider an equilibrium state. In this connection, we remark that 
the optimal trajectory provides a dynamical alternative to the usual prescrip- 
tion of equilibrium thermodynamics to calculate the free energy variation via 
a quasi static transformation. 

In Section [3] we discuss the notion of equilibrium from the standpoint of 
nonequilibrium under the assumptions formulated is Section [21 In particular 
we introduce the notion of macroscopic reversibility and discuss its relation- 
ship, actually equivalence, to thermodynamic equilibrium. An important 
outcome of our analysis is that the correct definition of equilibrium, for the 
driven diffusive systems here considered, is the vanishing of the currents. 
Nonetheless, in presence of external (electric) fields and boundary reservoirs 
an equilibrium state can be highly inhomogeneous. An example of such a 
situation is provided by sedimentation equilibrium in gravitational and cen- 
trifugal fields. In spite of this, the free energy is a local function of the state 
variables and coincides locally with the equilibrium free energy in absence of 
external fields and boundary driving. In particular there are no macroscopic 
correlations. 

In Section H] we derive from the Hamilton- Jacobi equation the generic 
existence of long range correlations in nonequilibrium states and we dis- 
cuss the conditions on the transport coefficients for their appearance in the 
two point correlation function. We also establish a simple criterion to de- 
termine whether density fluctuations are positively or negatively correlated. 
Finally, from the Hamilton- Jacobi equation, by expanding around the sta- 
tionary state, we derive a recursive equation for the correlation functions of 
any order. 

In Appendix |X] we discuss a specific example, the ABC model [10, 14] 
which does not satisfy all the assumptions formulated in Section [21 In fact 
the ABC model has the peculiarity to be a reversible system with a nonlocal 
free energy. We show nonetheless that this free energy can be easily derived 
from the Hamilton- Jacobi equation. This model has been brought to our 
attention by D. Mukamel and J.L. Lebowitz. 



2 Macroscopic systems out of equilibrium 

We here introduce the thermodynaniic description for out of equihbrium 
driven diffusive systems which are characterized by conservation laws. For 
simphcity of notation, we restrict to the case of a single conservation law, 
e.g. the conservation of the mass. The system is in contact with boundary 
reservoirs, characterized by their chemical potential Aq, and under the action 
of an external field E. We denote by A C R'^ the region occupied by the 
system, by x the macroscopic space coordinates and by t the time. We next 
state our basic axioms and their main implications. 

1. The macroscopic state is completely described by the local density p = 
p{t, x) and the associated current j = j{t, x). 

2. The macroscopic evolution is given by the continuity equation 

dtp + \/-j = (2.1) 

together with the constitutive equation 

j = J{p) = -D{p)Vp + x{p)E (2.2) 

where the diffusion coefficient D{p) and the mobility x{p) ^^^ d x d 
positive matrices. The transport coefficients D and x satisfy the local 
Einstein relation 

D{p)=x{p)f'M (2.3) 

where /o is the equilibrium free energy of the homogeneous system. The 
equations (12.11) - (12.21) have to be supplemented by the appropriate bound- 
ary conditions on dA due to the interaction with the external reservoirs. 
Recalling that Xo{x), x G dA, is the chemical potential of the external 
reservoirs, these boundary conditions are 

/^(p(x))=Ao(x) xedA (2.4) 

We denote by p = p{x), x E A, the stationary solution, assumed to be 
unique, of ^7Q), (Q, and fl^ . 

In postulating the constitutive equation (12. 2p we assumed a linear re- 
sponse with respect to the external field E and the thermodynamic forces. 
On the other hand the transport coefficients D and x ^i^re allowed to depend 
on the density p. In the case the system has more than one component, 
say n, the diffusion coefficient D and the mobility become nd x nd matrices. 
Moreover, in view of Onsager reciprocity, the matrix x is symmetric both in 



the space and in the component indices while D is symmetric only in the 
space indices. In such a case the local Einstein relation fl2.3p is D = xR 
where Rij = dp-dp. f^ does not depend on the space indices. In the context of 
stochastic lattice gases, (12. ip and (12. 2p describe the evolution of the empirical 
density in the diffusive scaling limit, see e.g. [6,18,26]; the validity of the local 
Einstein relationship (12.30 can be deduced from the local detailed balance of 
the underlying microscopic dynamics, see e.g. [26]. However, as we will see 
later, (12.30 it is also a consequence of the locality of the free energy and our 
third postulate once equilibrium is characterized within our theory. 

To state the third postulate, we need some preliminaries. Consider a time 
dependent variation F = F{t, x) of the external field so that the total applied 
field is E + F. The local current then becomes j = J^{p) = J{p) + x{p)F. 
Given a time interval [Ti, T2], we compute the energy necessary to create the 
extra current J^ — J and drive the system along the corresponding trajectory: 

L[T,,T.m= rdt{[j^ip^)-Jip^)].F)= rdt{F.xip^)F) (2.5) 

JTi JTi 

where ■ is the scalar product in R"', (■) is the integration over A, and p^ is 
the solution of the continuity equation with current j = J^{p). 

We define a cost functional on the set of space time trajectories as follows. 
Given a trajectory p = p(t, x), t G [Ti, T2], x G A, we set 

hTi,T2]ip) = 7 ^ in/ , kT^,T2]{F) (2.6) 

^ F : p'' =p 

namely we minimize over the variations F of the applied field which produce 
the trajectory p. The introduction of this functional in the context of fluid 
dynamics equations appears to be new: it is not the total energy provided by 
the driving field F but only the energy necessary to create the extra current. 
In fact, / (with the factor 1/4) has a precise statistical interpretation within 
the context of stochastic lattice gases: it gives the asymptotics, as the number 
of degrees of freedom diverges, of the probability of observing a space time 
fluctuation of the empirical density [1,2,5,18,19]. Note indeed that if p 
solves the hydrodynamic equation (12.10 . (12. 2p . (12. 4p its cost vanishes. We 
shall explain below the choice of the factor 1/4 in (12.60 . 

We next show that the optimal F in (12. 6p is given by F = VII where 
11: [Ti, T2] X A ^- M is the unique solution to the Poisson equation 

-v-[x(p)vn] =dtp + v-j{p) 

which vanishes at the boundary of A for any t G [Ti, T2]. Indeed, by writing 
F = Vn + F for each t G [Ti, T2] we have the orthogonal decomposition 

{F ■ x{p)F) = (vn ■ x(/3) vn) + {f ■ x{p)f) 



where we used that dtp + V ■ J^ (p) = 0. The above equation clearly shows 
that F = is the optimal choice. Hence 



i[T.,nip) = i X '^K [^^P + ^ ■ ^(z^)] ^^py' [^^P + "^ • ^('^)] ) (2.7) 

where the positive operator K{p) is defined on functions -u : A — i> M vanishing 
at the boundary dA by K{p)u = —V ■ (x(p)Vm). 

Our third postulate then characterizes the free energy J-'{p) of the system 
with a density profile p = p{x), a; G A, as the minimal cost to reach, starting 
from the stationary profile p, the density profile p, in an infinitely long time 
interval. 

3. The nonequilibrium free energy of the system is 

^(p)= inf /[_oo,o](p) (2.8) 

p:p(-oo)=p 

m=p 

As we shall see later, this is in fact proportional to the total work done 
by the optimal external field F along the optimal time evolution to reach the 
density profile p. The functional T has a precise statistical interpretation 
within the context of (non equilibrium) stationary states of stochastic lattice 
gases: it gives the asymptotics, as the number of degrees of freedom diverges, 
of the probability of observing a static fiuctuation of the empirical density [1, 
2,5]. In this sense the functional JF is the correct extension of the equilibrium 
free energy to non equilibrium thermodynamics. 

As shown in [1,2] the functional JF is the maximal solution of the infinite 
dimensional Hamilton- Jacobi equation 

where, for p that satisfies (12. 4p . bT jbp vanishes at the boundary of A. At the 
macroscopic level this condition refiects the fact that we consider variations 
of the density that do not change the boundary values. In the context of 
stochastic lattice gases, the boundary condition (12.41) is realized by superim- 
posing to the bulk dynamics suitable birth and death processes at the bound- 
ary; the probability of deviations of the empirical density at the boundary 
is then shown to be super-exponentially small in the scaling limit [2]. Note 
however that in the stationary state the density is not fixed at the boundary, 
i.e. ^(p) < cxD even if the profile p = p{x) does not satisfy the boundary 
condition ([231). 
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The arbitrary additive constant on the maximal solution of (12.91) is de- 
termined by the condition J^{p) = 0. By maximal solution we mean that any 
solution to (12.91) (satisfying F{p) = 0) is a lower bound for JF, as defined in 
(12.81) . This bound is sharp. Indeed, by considering the functional in (12.7J1 
as an action functional in variables p and dtp and performing a Legendre 
transform, the associated Hamiltonian is 

n{p, n) = (vn ■ x(p)vn\ + (vn • j(p)\ (2.10) 



where the momentum U vanishes at the boundary of A. By noticing that 
the stationary solution of the hydrodynamic equation corresponds to the 
equilibrium point (p, 0) of the system with Hamiltonian Ti and 7i(p, 0) = 0, 
the equation (12.91) is the Hamilton- Jacobi equation T-C[p,5J-'/5p) = 0. 

The optimal trajectory p* for the variational problem (12.81) is character- 
ized as follows. Let 

J*(p) = -2x(p)V^ - J(p) (2.11) 

dp 

then p* is the time reversal of the solution to 

dtp + v- r{p) = dtp + v- [d{p)Vp - x{p) [e + 2 v— ] } = (2.12) 

with the boundary condition (12.41) . We refer to [1,2,5] for the microscopic 
interpretation of J* in terms of the time reversed dynamics. In particular, 
the optimal external field F forcing the system on the time evolution p*(t) 
is so that J^{p) = —J*{p), namely 

F = 2V^ (2.13) 

dp 

This characterization of p* is obtained as follows. Let T be the maximal 
solution of the Hamilton- Jacobi equation and J* as in (12.111) . Fix a time 
interval [Ti,T2] and a path p(t), t G \Ti,T2\. We claim that 

%,T.](p)=^(p(T2))-^(p(Ti)) 

as can be shown by a direct computation using (12.71) . the Hamilton- Jacobi 
equation (12. 9p and the definition (12. lip of J*. From the identity (I2.14p we 
immediately deduce that the optimal path for the variational problem (12. 8p 
is the time reversal of the solution to (I2.12p . 



The classical thermodynamic setting considers only spatially homoge- 
neous system in which p does not depend on x. In this case the variation 
of the free energy between p and p is the minimal work required to drive 
the system from p to p, which is realized by a quasi static transformation, 
see e.g. [20]. For an equilibrium inhomogeneous density profile p, by a quasi 
static transformation we understand the following: pick a macroscopic point 
X and consider a macroscopically small volume around x, perform quasi static 
transformation from p{x) to p(a;), repeat independently for all x and sum the 
results. Of course, if the system is in contact with a reservoir also the con- 
tribution describing the interaction with the reservoir has to be included. 

We next show that according to definition (12.81) . the free energy T{p) is 
proportional to the work done by the optimal external field F on the system 
along the optimal trajectory p* . We emphasize that this is not a quasi static 
transformation in the above sense, but it is the solution of the hydrodynamic 
equation perturbed by the optimal external field. As we shall discuss later, 
the definition (12.81) of the free energy (with the factor 1/4 in (12. 6p ) agrees 
with the thermodynamic one for equilibrium states, but in general, due to 
presence of long range correlations, the variation obtained along a quasi static 
transformation as defined above gives a different result. 

Indeed, since p* is the time reversal of the solution to (12.121) . using (I2.13P 
we deduce 

t{p)-Hp) = / dt{—,dtp 

"0 



Sp 
£,,(|,V.J-(p-,)4/\*((-J.(p.),.F 



Namely, the variation J-'{p) — J^{p) is proportional the work done by the 
external field F along the optimal trajectory p* . The factor 1/2 above is 
connected to our definition of the transport coefficients and the associated 
Einstein relation (12. 3p . 

As shown in the next section, equilibrium states are characterized by 
J = J* . In this case the optimal path for the variational problem (12. 8p is the 
time reversal of the solution to the original hydrodynamic equation (12. ip - 
(12. 2p . In [2] we called such a symmetry property the Onsager-Machlup time 
reversal symmetry. 

3 Characterizations of equilibrium states 

We define the system to be in equilibrium if and only if the current in the 
stationary profile p vanishes, i.e. J{p) = 0. A particular case is that of a 
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homogeneous equilibrium state, obtained by setting the external field E = 
and choosing a constant chemical potential potential at the boundary, i.e. 
Ao(a^) = A. Let p = const, be the equilibrium density, i.e. p solves A = /o(p). 
It is then readily seen that the functional JF defined in (12. 8p is given by 



^(p) = / dx {fo{p{x)) - /o(p) -X[p{x) - p] } 



(3.i: 



in which the first difference is the variation of the free energy /o while the 
second term is due to the interaction with the reservoirs. The variational 
definition (12. 8p of the functional JF gives, when applied to equilibrium states, 
the variation of the equilibrium free energy (13. ip . This shows that the choice 
of the factor 1/4 in (12. 6p is appropriate. 

We next show that also for a non homogeneous equilibrium, characterized 
by a non constant stationary profile p(x) such that J{p) = the free energy 
functional JF can be explicitly computed. Let 



fip,x) 



dr / dr' ar') = U{p) - /o(p(x)) - f,{p{x)) [p - p(x)] 

p{x) J p{x) 



we claim that the maximal solution of the Hamilton- Jacobi equation is 

J'(p) = / dx f{p{x),x) 



Indeed 



5J^ 

5p{x) 



/o(p(^))-/o(p(^)) 



(3.2) 



(3.3) 



so that, by an integration by parts, 

V[/^(p)-/^(p)]-x(p)V[/^(p)-/^(p-)] 

+ {[Up)-fM]^-[D{p)Vp-x{p)E 

'y[ap)-ap)\-x{p)\^ap)-E 







where we used (12. 3 p and Vfoip) — E = —x{p)^^J{p) = 0. Therefore the 
functional JF in (13.20 satisfies the Hamilton- Jacobi equation (12. 9p . To show 
it is the maximal solution, recalling (12.70 . simple computations show that 



'[Ti,T2 



(P) 



^(/3(T2))-^(p(Ti)) 



+ 



IJ'dt^ [dtp - V ■ Jip)] Kip)-' [dtp - V ■ J{p)] 



(3.4) 



which clearly implies the maximality of JF. 

We emphasize that the above argument depends crucially on the structure 
(12.21) of the current and on validity of the local Einstein relation fl2.3p . In 
Appendix |X] we discuss the example of the ABC model [10, 14] in which 
these conditions are not met and the free energy is not a local functional 
even if J[p) = 0. Another example is given by the anisotropic zero-range 
process [9,27] for which - as discussed in [27] - the local Einstein relation is 
violated. 

We next show that the condition J{p) = is equivalent to either one of 
the following statements. 

- There exists a function A : A ^ M such that 

E{x) = VA(x) , X G A A(x) = Ao(x) , x e dA (3.5) 

- The system is macro scopically reversible in the sense that for each 
profile p we have J*{p) = J{p)- 

We emphasize that the notion of macroscopic reversibility does not imply 
that an underlying microscopic model satisfies the detailed balance condition. 
Indeed, as it has been shown by explicit examples [15,16], there are non 
reversible microscopic models which are macroscopically reversible. 

We start by showing that J{p) = if and only if (13. 5p holds. From the 
local Einstein relation (12. 3p and J{p) = we deduce 

E{x) = f;;{p{x))vp{x) = vap{x)) 

hence (13. 5p . Conversely, let the external field E be such that (13. 5p holds. 
Since /q' is positive the function /q is invertible and we can define p{x) = 
(/q)^^(A(x)). The profile p satisfies (12. 4p as well as J{p) = 0, consequently 
also V ■ J{p) = so that it is the unique stationary profile. 

We next show that J{p) = if and only if J(p) = J*{p)- By evaluating 
the Hamilton- Jacobi equation for p = p we deduce V[5J-'{p)/5p'\ = 0. From 
(12. lip we then get J{p) = 0. To show the converse implication, note that if 
J{p) = then JF is given by (13.20 . We deduce 

x(p)V^ = Dip)Vp - x{p)E = -J{p) 
dp 

hence, recalling (E^D and (EH]), J{p) = J*{p). 

We also note that macroscopic reversibility J{p) = J*{p) implies the in- 
variance of the Hamiltonian 7i in (I2.10p under the time reversal symmetry, 
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see [2] , (p, n) \-^ (p, 6J-'/6p — U) , where J-' is the maximal solution of the 
Hamilton- Jacobi equation fl2.9l) . As remarked before, macroscopic reversibil- 
ity implies the Onsager-Machlup time reversal symmetry, on the other hand 
we shall give at the end of this section an example of a nonequilibrium sys- 
tem for which such a symmetry holds. Hence macroscopic reversibility is a 
stronger condition. 

So far we have assumed the local Einstein relation and we have shown 
that -for equilibrium systems- it implies (13.21) . Conversely, we now show that 
macroscopic reversibility and (13.21) imply the local Einstein relation (12. 3p . By 
writing explicitly J{p) = J*{p) we obtain 

- [xip)Rip) - D{p)]Vp + xip) [Rip) - X-\p)D{p)\V-p = (3.6) 

where R is the second derivative of /o in the case of one-component systems 
while Rij = dp-dp^fo for multi-component systems. In (13. 6p we used, be- 
sides (13.21) . J{p) = to eliminate E. Note that J{p) = follows from the 
Hamilton- Jacobi equation and J{p) = J*{p) without further assumptions. 
Since p and Vp are arbitrary the local Einstein relation D = xR follows 
from (13.61) . 

We have defined the macroscopic reversibility as the identity between the 
currents J{p) and J*{p)- We emphasize that this is not equivalent to the 
identity between V ■ J{p) and V ■ J*{p)- Indeed, we next give an example of 
a non reversible system, i.e. with J{p) ^ 0, such that the optimal trajectory 
for the variational problem (12.81) is the time reversal of the solution to the 
hydrodynamic equation (12. II) .( 12121) . and (12. 4p . 

Let A = [0,1], -D(p) = x(p) = 1, Ao(0) = Ao(l) = A, and a constant 
external field -E 7^ 0. In this case hydrodynamic evolution of the density 
is given by the heat equation independently of the field E. The stationary 
profile is p = A, the associated current is J{])) = E =/^ Q. By a computation 
analogous to the one leading to (13.41) . we easily get that 

and the optimal trajectory for the variational problem (12. 8p is the time rever- 
sal of the solution to the heat equation. On the other hand J{p) = — Vp + E 
while J*{p) = -Vp- E 

We remark that, even if the free energy JF is a non local functional, the 
equality J(p) = J*{p) implies that the thermodynamic force V6J-'/6p is local. 
Moreover, the Hamilton- Jacobi equation reduces to the statement 

J{p) = -X(P)V^(P) (3.7) 
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This identity, viewed as a special case of fl2.11l) . represents the general form 
- for reversible systems - of the relationship between currents and thermo- 
dynamic forces. It holds both when the free energy is local and nonlocal. 
For reversible systems it is therefore enough to integrate (13. 7p to compute 
the free energy. An example of such a situation with a nonlocal free energy, 
illustrated in Appendix \^ is provided by the ABC model on a ring with 
equal densities [10,14]. 

4 Correlation functions 

In classical equilibrium thermodynamics the second derivatives of the the 
thermodynamic potentials like the free energy are response coefficients like 
the compressibility. In a statistical mechanics framework they are identified 
with the Fourier transform evaluated at zero of correlation functions of the 
underlying Gibbs measure. Analogous interpretation is given to the higher 
order derivatives. 

As we defined the (nonequilibrium) free energy as a functional on the 
set of density profiles, we can obtain the (nonequilibrium) density correla- 
tions functions of arbitrary order in terms of the functional derivatives of 
JF. In general the functional JF cannot be written in a closed form, but - 
by a suitable perturbation theory on the Hamilton- Jacobi equation (12. 9p - 
we can derive such correlations functions. In this section we first discuss the 
two-point correlation and establish a criterion to decide whether the den- 
sity fiuctuations are positively or negatively correlated. Then, we deduce a 
recursive equation for the correlation function of any order. 

We remark that we are concerned only with macroscopic correlations 
which are a generic feature of nonequilibrium models. Microscopic correla- 
tions which decay as a summable power law disappear at the macroscopic 
level. 

We introduce the pressure functional as the Legendre transform of free 

energy JF 

gih) = snp{{hp)-J^ip)} 
p 

By Legendre duality we have the change of variable formulae h = -^, p = -^, 
so that the Hamilton- Jacobi equation (12. 9p can then be rewritten in terms 
of Q as 

>-x(f)v.)-(v...(f)vf-.(|).)^0 ,4.) 

where h vanishes at the boundary of A. As for equilibrium systems, Q is the 
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generating functional of the correlation functions. In particular by defining 

Cix V) = ^'^(^) 
^ '^^ 6hix)6h{y) h=o 

we have, since J-" has a minimum at p, 

1 



g{h) = {h,p) + ^{h,Ch) + o{e) 



or equivalently 



Hp) = l{ip-p), C-\p - -p)) + o{{p - -pf) 

By expanding the Hamilton- Jacobi equation fl4.ll) to the second order in 
/i, and using that 6Q/6h{x) = p{x) + Ch{x) + o{h'^), we get the following 
equation for C 

Vh ■ [x{p)Vh - V{D{p)Ch) + ^{p){Ch)E\ ) = (4.2) 

We now make the change of variable 

C(x, y) = Ceq(a;)5(a; - y) + B{x, y) 

where Ceq(x) is the equilibrium covariance. By using (12. 3p we deduce that 

aq(x) = D-\-p{x))x{-p{x)) 

Equation (14.21) for the correlation function then gives the following equation 
for 5 

C)B{x,y)=a{x)5{x-y) (4.3) 

where £^ is the formal adjoint of the elliptic operator C = L^ + Ly given by, 
using the usual convention that repeated indices are summed, 

L. = D,,ipi^))d.A, + ^,^{-p{x))E,{x)d.^ (4.4) 

and 

«(^) = d.. [^^M^)) D-^{-p{x))Ux)\ 

where we recall J = J{p) = —D{p{x))'Vp{x)+x{p{x))E{x) is the macroscopic 
current in the stationary profile. 

In agreement with the discussion in Section [HI for equilibrium systems 
J = so that we have a = 0, hence 5 = 0, namely there are no long range 
correlations and C{x,y) = Ceq{x)6{x — y). Moreover, since C is an elliptic 
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operator (i.e. it has a negative kernel), the sign of B is determined by the 
sign of a: if a{x) > 0, Vx, then B{x,y) < 0, Wx,y, while if a{x) < 0, Vx, 
then B{x,y) > 0, Wx,y. For example, consider the following special case. 
The system is one-dimensional, d = 1, the diffusion coefficient is constant, 
i.e. D{p) = Do, the mobility x{p) is a quadratic function of p, and there is 
no external field, E = 0. Then 

B{x,y) = --^x"{^p)'^-\x,y) (4.5) 

ZJJq 

where A~^{x,y) is the Green function of the Dirichlet Laplacian. Two re- 
markable models, the symmetric exclusion process, where x{p) = p(l ^ p); 
and the KMP process, where x{p) = P^? meet the above conditions. Then 
04.51) shows that their correlations have opposite signs. 

We next derive, using the Hamilton- Jacobi equation 04.11) . a recursive 
formula for the ra-point correlation function Cn{xi, . . . , x„). This is defined 
in terms of the pressure functional Q as 

Cn{Xi, ...,Xn)= ^,, . JT-f 7 (4-6) 

0h[Xi) ■ ■ ■ 0h[Xn) h=0 

so that Ci{x) = p{x) and C2 is the two-point correlation function discussed 
above. Note that for equilibrium systems the free energy is given by 03.21) . 
hence we have no long range correlations and the n-point correlation function 
is 

Cn,eqiXi, . . . , X„) = C„(p(Xi) )(5(Xi - X2) " " • (5(Xi - X„) (4.7) 

where c„(p) = £L(/^)-i(x)|^^^,(_^. 

By expanding the functional derivative of Q we get 

||)=*.) + E^ft.(/..0 (4.8) 

where 

Qn{h;Xi)= dX2...dXn+lh{x2) ...h{Xn+l)Cn+l{Xi,X2,...,Xn+l) (4.9) 

Ja 
If we consider the terms of order /i"+^ in the Hamilton- Jacobi equation 04. ip 
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we get, using 



V/i- 



E 



and after some algebraic manipulations 
1 



T,k'^ik=1T—l 



ii!(l!)^ii2!(2!)^^ 



-X 



iT.k->-k)(n 



{p){Qir{Q2r...vh 



- E 

*l,«2,--->0 

+ E 

il,42,--->0 



tAi\yH2\{2\\ 



«2 . . . V 



vfz}(E. •'.-!) (p)(^^)n(^^ 



ii\{l\YH2\{2\Y^ 



-X 



(T,kik)(n 



{p){Gir{<52r---E 



where, for / given by either x or -D, the notation f^^'' stands for the /c-th 
derivative of /. 

Since the above identity holds for every function h vanishing at the bound- 
ary of A, we get from it a recursive formula for the n-point correlation func- 
tions which generalizes fl4.3p . We start by introducing some notation. For 



Xkl 



n > 1, we let Cn^i be the formal adjoint of the operator £„+i = X]fc=i -^ 
where L^ is defined in (14.41) . Given functions / and g in n + 1 and m + 1 
variables respectively, their 1st shuffle product f^ig is, by definition, the 
following function in m + n + 1 variables 



/#l5'(3^1) • • • 5 Xm+n+l) 

= 7 — \ — ^Ty]/(^i' 

[m + ny. ^^-^ 



• Xn„ + i)9[^lj ^TTn 



•^Tm + n + lJ 



where the sum is over all permutations of the indices {2, 3, . . . , m -|- n + 1}. 
Clearly #i is a commutative associative product, and fi^ig is symmetric in 
the variables X2, ■ ■ ■ , Xm+n+2- Furthermore, given a function f in n variables, 
we denote by f^^"^ its symmetrization, namely 

/ [Xl, . . . , Xn) = I / J jy^TTl ) • • • ; ^TTn) 



TreSn 



Given a sequence X= {ii,i2, ia, ■ ■ ■) oi non negative integers with only finitely 
many non zero entries, we also denote 



E(.-) 
N{i) 



ii + i2 + ■ 
k + 2^2 + 



E 



Ik 



^kiu 



K{l) = i,\{l\rH2\{2\) 



W^kW^ 



*fc 
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Finally, for Twitli NlV) = n, we let 

C^x,, . . . ,x„+i) = (Cf '^#iC3#^^^#i . ■ ■)(xi,X2, . . . ,a;„+i) 

Then the recursive equation satisfied by the n + 1-point correlation function 
reads 



[n + l 



-—-C^^^-^Cn+liXi,X2, ..., Xn+l) (4.10) 



{ Yl Wt)^'-- ■ {x^'^^^\p{Xl))CA^Xi, . . . ,X„)V,.,5(Xi - Xn+l) 



N{i)=n-1 



E ;^V., ■ V.,(D(^«-i)(p(xi))Cr(xi,...,x„+0) 



N{L)=n,i„=0 



sym 



Af(r)=n,i„=0 

It is easy to check that, for ra = 1, equation fl4.10p reduces to (14.31) . Moreover, 
in the case of equilibrium states the solution to (I4.10p is (14.71) . We refer 
to [2, § 4.5] for the application to the boundary driven simple exclusion 
process. 



A The ABC model 

We here consider - both from a microscopic and macroscopic point of view 
- a model with two conservation laws. Given an integer A^ > 1 let Z^r = 
{1, . . . , N} be the discrete ring with A^ sites so that A^ + 1 = 1. The micro- 
scopic space state is given by VLn = {A, B, C}^^ so that at each site x G Z^v 
the occupation variable, denoted by rj^, take values in the set {A, B, C}; one 
may think that A, B stand for two different species of particles and C for 
an empty site. Note that this state space takes into account an exclusion 
condition: at each site there is at most one species of particles. 

We first consider a weakly asymmetric dynamics that fits in the frame- 
work discussed in Section [2] that is defined by choosing the following transi- 
tion rates. If the occupation variables across the bond {x,x + 1} are (^,C)) 
they are exchanged to (C, with rate c^x_^_^ = exp{{E^ — E(^) / {2N)} for fixed 
constant external fields Ea,Eb,Ec- This choice satisfies the so-called local 
detailed balance condition that, in general, requires the weakly asymmetric 
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rates c^x_^i to satisfy 

<.+i(r/^'^+^) = cf;,.+i(r/) (A.l) 

X exp I V^,^+i H{ri) - tt^ {Eax1a{Vx) + Ebx1b{Vx) + EcxldVx)) 

X 

where r^^'^'+i is the configuration obtained from rj by exchanging the occupa- 
tion variables in x and x + 1, 'Wx,x+if{fl) = f(^rj^'^~^'^) — f[ri), and H{ri) is the 
energy of the configuration rj, which is constant in this modeL 

Following [12,18,19], the hydrodynamic equations for the densities of A 
and B particles are given by 

d ( ^^\ = a( ^^\ -V ^ ^^^^ ~ ^^^ -PAPB \ ( Ea-Ec 



Pb J \ Pb ) V ~PaPb Pb(1 - Pb) J \ Eb - Ec 

(A.2) 
of course the density of C particles is then pc = I — Pa — Pb- 
The functional /[ri,T2] in (12.71) with D = 1 and mobility 

V -PaPb Pb{^ - Pb) J 

is the dynamical large functional associated to this model. The free energy 
is the maximal solution of the Hamilton- Jacobi equation (12. 9p which can be 
easily computed. Namely, 



^mA,mB {PA,Pb)= I dx PA log ^ + pB log 
J I niA 



Pa , ^ ^^^PB_ 

rriB 



+ {1-PA- pb) log / ^^ ^^ 1 (A.4) 

1 — TTiA — rns-i 

where Jdx pa = rriA and Jdx pb = trb- If Ea, Eb and Ec are not all equal, 
this model is a nonequilibrium model nevertheless, in view of the periodic 
boundary conditions, its free energy is independent of the external field; 

see [6, § 3.4]. 

We next discuss a different choice of the weakly asymmetric perturbation 
which, as we shall see, does not fit in the scheme discussed in Section [21 This 
choice is the one referred to in the literature [10,14] as the ABC model. The 
transition rates are the following. If the occupation variables across the bond 
{x, X + 1} are (^,C); they are exchanged to (CO with rate exp{V^(C''7)/iV} 
where V{A,B) = V{B,C) = V{C,A) = -13/2 and V{B,A) = V{C,B) = 
V{A, C) = (3/2 for some /3 > 0. Therefore the A-particles prefer to jump 
to the left of the 5-particles but to the right of the C-particles while the 
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i?-particles prefer to jump to the the left of the C-particles, i.e. the preferred 
sequence is ABC and its cydic permutations. These rates do not satisfy the 
local detailed balance ( lA.ip . 

Again by the methods developed in [12,18,19], the hydrodynamic equa- 
tions are 

\ PB J V Mpa,Pb) J 

where 

T/^ ^^ ( Ja{pa,Pb)\ ( -Vpa + I3pa{1-'^Pb- Pa)\ ^ . r.-. 
Apa, Pb) - ^ j^^p^^ Pb))-\ -Vps + f3pB{2pA + PB-1) ) ^^-^^ 

The asymmetric term in the hydrodynamic equation ( 1A.5I1 is not of the form 
V ■ (x(p)-E') as in f l2.1l) -( l231) . Hence the results obtained in Section [3] do not 
apply. 

The appropriate cost functional is however still given by (12.71) with J as 
in (1A.6P and x cls in (lA.Sp . The solution of the Hamilton- Jacobi equation 



(12.91) then gives the free energy. In the case of equal densities j dx pA 
Jdx Pb = 1/3, a straightforward computation shows that for any positive /3 
the solution is given by the functional 

^1 1 {pa, Pb) =^11 (pa, Pb) 



+(3 dx dyy^pA{x)pB{x + y)+pB{x)[l-pA{x + y)-pB{x + y)] 

+ [1 - Pa{x) ~ pb{x)]pa{x + y)^ + X (A.7) 

where J-'l i is the functional in flA.4p with ttia = ttib = 1/3 and x is the 

3'3 

appropriate normalization constant. This result has been already obtained in 
[10] by direct computations from the invariant measure. Indeed, in this case, 
the ABC model is microscopically reversible and the invariant measure can 
be computed explicitly. From a macroscopic point of view, the reversibility 
of the model is expressed as the identity, which holds in the case of equal 
densities, J{p) = J*{p). Hence, as discussed at the end of Section [31 the free 
energy (JA.71) can also be macroscopically derived by integrating (13.71) . 

We mention the recent preprint [8] , of which we became aware after the 
completion of this work, in which the nonequilibrium free energy of the ABC 
model is also computed from the variational principle (12.81) . 
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